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ABSTRACT 

A   simple    two- level   numerical  model   using    the   quasi-geostrophic 
forecast   equations    is   developed.      Equations    are    linearized   and 
friction   is    introduced    in   the    surface    layer.      Solutions   are    obtained 

numerically  by   using    the    initial   value   approach.      Two  wind    profiles, 

2 

U  =  -U  tanh  y/y  and  U  =  U  sech   y/y  ,  are  used  and  these  are  known 

o      J    J o  o         J  o 

to  be  unstable.   For  each  wind  profile  the  growth  rate  is  determined 
as  a  function  of  the  wave  number.   Some  observed  features  of  easterly 
waves  are  reproduced  in  the  numerical  solutions. 
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1 .   INTRODUCTION 

Synoptic  scale  disturbances  in  the  tropical  atmosphere  have  long 
been  recognized  (Palmer,  1951;  Riehl,  1954).   In  general  these  dis- 
turbances possess  a  "cold-core"  structure  with  the  colder  air  and  most 
active  weather  located  east  of  the  wave  (Yanai  and  Nitta,  1968b). 
Vertical  motions  associated  with  these  waves  are  generally  upward  east 
of  the  wave  axis  and  downward  to  the  west  (Yanai  and  Nitta,  1967).   The 
source  of  energy  for  these  waves  is  not  fully  known,  however,  Charney 
(1963)  argues  that  when  condensation  is  absent  the  large  scale  motions 
tend  to  be  quasi-horizontal  and  quasi-nondivergent .   If  this  argument 
is  correct  then  the  barotropic  instability  that  arises  as  the  result 
of  horizontal  shear  in  the  Intertropical  Convergence  Zone  could  be  the 
source  of  energy  for  the  disturbances  (Schminke,  1968). 

Kuo  (1949)  studied  the  stability  characteristics  of  a  barotropic 
zonal  current,  and  showed  that  the  barotropic  zonal  current  is  stable 
if  the  gradient  of  absolute  vorticity  has  the  same  sign  in  the  zonal 
current.   Charney  (1964)  discussed  the  case  where  horizontal  shear 
instability  is  produced  in  the  equatorial  convergence  zone  when  the 
zone  is  located  away  from  the  equator  and  the  converging  air  masses 
carry  their  angular  momentum  with  them. 

The  purpose  of  this  paper  is  to  devise  and  test  a  numerical  model 
that  is  structured  realistically,  that  is,  the  perturbation  has  a  con- 
figuration and  behavior  similar  to  observed  easterly  waves.   The  main 
feature  desired  is  to  have  the  associated  weather,  and  therefore  the 
vertical  motion,  upwind  of  the  trough  axis. 


2.   BASIS  OF  THE  MODEL 

The  assumption  that  the  source  of  energy  for  the  wave  mechanism  is 
barotropic  instability  leads  to  the  selection  of  a  basic  current  that 
is  independent  of  height.   Baroclinicity  has  been  omitted  but  the 
atmosphere  remains  stratified.   Charney  (1963)  used  a  scale  analysis  to 
show  that  the  barotropic  vorticity  equation  governs  tropical  motions  in 
the  absence  of  condensation.   In  another  paper  Charney  (1969)  discusses 
the  role  of  the  phase  velocity  in  energy  propagation  in  an  easterly 
regime.   Following  Schminke  (1968)  the  assumption  is  made  that  the 
easterly  waves  can  be  roughly  described  by  a  two-level  quasi-geostrophic 
model.   It  is  recognized  that  the  quasi-geostrophic  approximation  is  not 
very  accurate  close  to  the  equator,  however,  the  zone  of  dynamic 
instability  is  assumed  to  be  associated  with  the  Intertropical  Conver- 
gence Zone  which  is  some  distance  from  the  equator  so  that  the 
approximation  is  relatively  accurate. 

The  type  of  wind  profile  chosen  will  play  a  large  role  in  determining 
the  stability  of  the  disturbance.   Jacobs  and  Wiin-Nielsen  (1966)  inves- 
tigated barotropic  instability  and  found  that  there  are  several  unstable 
modes.   Since  a  basic  zonal  current  with  no  change  in  the  vertical  is 
used  in  this  model,  the  wave  with  no  vertical  motion  would  predominate 
(Jacobs  and  Wiin-Nielsen,  1966).   Easterly  waves  are  observed  to  have 
vertical  motion,  so  a  method  must  be  found  to  introduce  it.   This  is 
accomplished  by  introducing  friction  in  the  surface  boundary  layer. 
Thus  vertical  motion,  which  depends  upon  the  stratification  of  the 
atmosphere,  is  now  forced  into  the  model. 


3.   THE  FORECAST  EQUATIONS 

A   simple    two-level  model   is    constructed   by  dividing   the   entire 
atmosphere    into    four    layers    of   constant    pressure   differential,    Ap/2 
(fig.    1),    numbered    0   to   4    from  top   to  bottom.      Vertical  motion   is 
assumed    to  be    zero   at    the 

0  ^^ 0 
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1 
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3 


2Ap 


00  =    -    pgw 


Fig.    1.      Two-level   model   used    for    prediction. 

top   of   the   atmosphere,   while    the   vertical   motion   term  at    the   earth's 

surface    is 

m  =  -   pgw4.  (3.1) 

Charney  and   Eliassen    (1949)    used    the   Ekman   theory   to  derive   an 

expression    for  w,    which    is  ,  i* 

w4  =j(^)  sin  2a  C4,  (3.2) 

o 

where  A   is  the  kinematic  eddy  viscosity  and  a,  the  inflow  angle.   The 
surface  geostrophic  vorticity  can  be  approximated  by 

C4  =  C3   =  V2Y3.  (3.3) 

Following    the    development    of   Thompson    (1961),    but    using   different 
notation,    begin  with   the    quasi-geostrophic   vorticity  equation 

d        2  i  2  rW  chi 

§-  V  Y  +     k  x  VY-V   (V  Y)  +  3     -F  -  f     TT  =  0.  (3.4) 

dt  '  o   ox  o   op 

Notice  that  when  oo  =  0,  the  above  equation  becomes  the  barotropic 
vorticity  equation,  which  is  sufficient  to  describe  the  main  instability 


Apply  this  equation  at  levels  1  and  3  giving 

n        2  2  ^1  ^2 

§-  V  Y,   +    |k  x  VY,'V(V  Y,)  +  0     -s—  -  f     -^  =  0,  (3.5) 

dt  1  1  1  o   dx  o  Ap 

d        2  ,  2  ^Y3  (<04    '  C°2) 

§T  V  Y,  +     k  x  V!    -V(V1J  +  0     vr+f  — "h ~  -  0.      (3.6) 

dt  3  3  3  o   ox  oAp 

Smaller    layers    could   be    used    if   it    is   desired    to  stop  at   some   intermediate 

point,    say   the    tropopause,    rather   than    including    the   entire   atmosphere. 

Next   consider   the   quasi-geostrophic    first    law  of    thermodynamics    in 

the    form 


Ie  &+'■"".  *•"(£>!-•-'».       <3-7> 


where 

2-  ■-<= 


•-VK+S-1  <"•> 


p  g 

When  this  equation  is  applied  at  level  2  the  result  is 

d  ,      (W  ApaD2 

ft"  (Y1  -  Y3)  +  Ik  x  V  -^  •  V(Yr¥3) -^   =  0.   (3.8) 

o 


Define  the  following  quantities 


Y 
m 

= 

*i 

+ 

Y 
T3 

2 

Y 

= 

*1 

- 

Y 
T3 

2 

(3.9) 
(3.10) 


which    implies    that 


Yl    =  Y™  +  Yt  <3-U> 

1  m  1 

Y3  =  Ym  "  V  (3'12) 

Here  YT  is  proportional  to  the  layer  thickness,  therefore  is  a  measure 

of  its  mean  temperature.   Using  these  definitions  add  (3.5)  and  (3.6) 

and  divide  the  result  by  2,  obtaining 

7\       ?  ■)  o  ^Y  m, 

%-  V^Y     +    |k  x  VY    -VOTY   )  +    |k  x  V^-V^V   )  +  P     TT~f     ^T-  =  0. 
dt  m        '  mm'1  T  o  dx  o  2^p 

(3.13) 
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For  the  second  forecast  equation  subtract  (3.6)  from  (3.5)  and  eliminate 
ecu  using  (3.8) 

|^  (V2-  [i2)    YT  +  |k  x  VYm-V(V2-^)  YT  +  |k  x  VYT. 


T 

oYP 


CD, 


2         W1T       4 
V(V^Y  )  +  p  rrr  +   f  oT"  =  o, 

m     o  ox     o  2Ap 


(3.14) 


where 


9   2f 
2     o 

M  = 


,  2 

Ap  a 


(3.14a) 


These  are  the  prediction  equations  for  the  model.   They  can  be 
linearized  by  separating  the  flow  into  an  east-west  current  that  varies 
only  in  y,  and  a  small  departure  from  this  flow.   In  this  event  it  is 
possible  to  treat  various  waves  in  x  independently,  hence  the  fields 
may  be  defined  as  follows 

Y  =  E(y)  +  A(y,t)  cos  kx  +  B(y,t)  sin  kx,       (3.15) 


YT  =  C(y,t)  cos  kx  +  D(y,t)  sin  kx, 


(3.16) 


where  k  is  the  x  wave  number. 


Substitute  expressions  for  Y  and  Ym  into  (3.13)  and  (3.14), 

m      T 

separate  the  various  sine  and  cosine  terms,  neglecting  all  products  of 
the  quantities  A  through  D,  and  equating  coefficients  of  the  cosine  kx 
terms  gives 


at 


^A 


Ak   =  k 


cfc   cHB 

^ay2 


f-2-^B 


3y" 


] 


B  Bk  -  K 
o 


X 


*' 


(A  -  C) 


For  the  sine  terms  the  result  is 


h<&-* 


=  k 


j*  Ak2  -  a  4  +  4  a 

Lay  By   oy2        oy3 


8      Ak 
o 


K  (Z-   -  k2 

V 


:)    (B    -   D). 


(3.17) 


(3.18) 
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Repeating    the    procedure    for    (3.16)    the    final   equations   are 


5£  ^—  -  ck 


dy 


,  z\     .   r0E  a  D      oE 


D    (k     +  |1   )    - 


a3E 


dy 


3   UJ 


\q  Dk  +  K  (-—  -  k2)   (A   -   C)  , 

dy 


(3.19) 


and 


V 


3t  V7T  "  Dk    "  D^  i  =  k 


oE    _    .,  2  2.         oE   ^2C 


Lay 


dy 


V 


^-|  C     +  eo  Ck  +  K  (^   -   k2)   (B    -   D)  , 


(3.20) 


dy 


V 


where 


K  = 


2RT        f  o 


1/2 


(3.21) 


It    is    now  necessary   to   derive    the   computational   equation   for   vertical 
motion   starting  with 


roY 


^2  ""  Apa 


^Yn 


^L 


T    ,  ~*T    ,  ~*T 

— <r h    U       — ^r h    V       — =r—     , 

at  m      ax  m      ayJ 


where 


u      = 
m 


,       aY 
i  m 


f        ay    ' 
o        J 

.      aY 
m 

m        f      ax 
o 


(3.22) 


(3.23) 


(3.24) 


Substituting  equations  (3.23),  (3.24),  (3.15)  and  (3.16)  into  (3.22), 
collecting  terms,  and  again  neglecting  all  products  of  the  quantities 
A   through   D   results    in 

2 


^2  "  Apa 


cos    kx  (j|  _  L.  4|  d)  + 

o 


Be 


■in  kx  (^  +  T  ^  ( 


(3.25) 


This  is  the  equation  used  in  the  model  for  the  computation  of  vertical 
motion.  The  prediction  equations  and  the  vertical  motion  equation  are 
now  in  the  form  to  be  used  in  the  model. 

12 


4.   BOUNDARY  CONDITIONS  AND  FINITE  DIFFERENCING  SCHEME 

The  finite  differencing  scheme  used  is  illustrated  below  with  a 
sample  variable  N; 


ri(Ni+l    "   Ni-l) 


a2N 

-  1  ( 

3y2 

V 

a3N 

1 

ay3 

2H3 

Ni+1  -  2Ni  +  Nl-l) 


Nl+2    "   2Vl  +  N 


(N.    -    2N.     ,    +  N.    .) 
l  l-l  1-2      J 


(4.1) 
(4.2) 


(4.3) 


where  i  is  the  grid  index  and  H  is  the  distance  between  grid  points. 

Centered  time  differences  are  used  for  all  quantities  except  those 
involving  friction.   The  frictional  terms  are  computed  at  time  (t  -  At). 
The  first  step  in  all  cases  is  a  forward  time  step.   The  second  order 
equations  for  time  tendencies  are  solved  by  the  exact  method  of  Richtmyer 
(1957,  p.  101). 

Yanai  and  Nitta  (1968a)  studied  the  problem  of  finite  difference 
approximations  in  solving  dynamic  instability  problems  of  non-divergent 
barotropic  currents.   They  showed  that  for  both  symmetric  and  anti- 
symmetric zonal  currents  the  exact  boundary  conditions  can  be  replaced 
by  rigid  boundary  conditions  at  a  distance  equal  to  the  half  width  of 
the  shearing  wind  belt.   At  any  boundary  distance  larger  than  this  value 
the  numerical  and  theoretical  values  tended  to  have  a  constant  difference 
Therefore  rigid  boundaries  are  placed  at  y  -  -w/2  and  y  =  w/2,  where  w  is 
the  total  width  over  which  computations  are  made.   With  the  exception  of 
E  (the  basic  current)  all  variables  have  the  boundary  condition  A  =  B  = 
C=D=0aty=  -w/2  and  at  y  =  w/2.   When  second  derivatives  are 
required  on  the  boundary  they  are  set  equal  to  zero. 
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Fig.    2.      Wind    profile   U   =    -U    tanh    y/y 

o  o 


5.   WIND  PROFILES  AND  INITIAL  CONDITIONS 

Schminke  (1968)  investigated  and  discussed  the  structure  of  easterly 

waves    from  a   similar   model   using   one  wind    profile    (U  =   -U,_tanh(y/y    )) 

U  o 

and   one    non-dimensional  wave    number    (l/2y   ).      It    is    proposed    to  expand 

2 
this  to  include  another  profile  (U  =  U  sech  y/y  ).   Additionally  the 

behavior  of  disturbances  at  various  non-dimensional  wave  numbers  under 

both  profiles  will  be  considered. 

Charney  (1963)  shows  by  scale  analysis  that  in  the  absence  of 
condensation,  flow  in  the  tropics  is  governed  by  the  barotropic 
vorticity  equation.   The  assumption  is  made  that  the  source  of  energy 
for  easterly  waves  must  be  barotropic  instability.   Since  potential 
energy  is  assumed  not  to  be  an  important  energy  source  the  basic 
temperature  field  may  be  set  equal  to  a  constant  and  is  invariant. 
However,  small  temperature  fluctuations  are  observed  with  easterly 
waves,  so  a  two  level  model  is  needed  (Schminke,  1968). 

If  the  waves  are  to  grow,  the  zonal  wind  profile  must  be  baro- 

tropically  unstable.   Following  the  suggestion  of  Charney  (1963)  that 

the  shear  layer  near  the  Intertropical  Convergence  Zone  may  be  formed 

by  transport  of  angular  momentum,  the  unstable  profile  chosen  to 

approximate  this  is 

U  -  -UQtanh  y/yQ,  (5.1) 

illustrated  in  figure  2.   This  seems  most  appropriate  during  monsoonal 

conditions.   The  second  profile  chosen, 

U  =  U  sech2y/y   ,  (5.2) 

o     J    J  o 

is  probably  more  applicable  to  mean  flow  conditions,  and  is  illustrated 
in  figure  3. 
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Garcia  (1956)  found  that  the  profile  given  in  (5.1)  is  barotropically 
unstable  when  0  <  k  <  y    .   Betchov  and  Criminale  (1967)  have  analyti- 
cally, and  by  use  of  the  eigenvalue  approach,  investigated  this  profile 
throughout  the  range  mentioned  above  and  found  the  maximum  instability 

at  a  non-dimensional  wave  number  of  0.45  (ky  =  0.45).   With  these 

J  o 

results    in  mind   this    paper    is   confined    to    the    range    0.25   to   0.60   for 

the    non-dimensional  wave    number  when  dealing  with    the   hyperbolic 

tangent    profile    (5.1).      Betchov   and   Criminale    (1967)    found   a    rather 

flat    zone    of  maximum   instability   for   the    hyperbolic   secant    squared 

profile    (5.2)    in   the   range    0.80   to    1.05    for    the    non-dimensional  wave 

number.      Some   computations   were   done    in   the    lower   range    (0.25   to   0.60), 

but   the  majority   of   the    experiments   with    this    profile  were    confined    to 

the    range    0.70   to    1.15.      Results    of  experiments   with   both   profiles   will 

be   discussed    in  a    later   section. 

With    the   wind    profile    independent    of  height   and  with    no    friction  Y 
r  m 

and    YT   should  behave    independently.      Jacobs    and   Wiin-Nielsen    (1966) 

showed   that   the   growth   rate    for   the    Y      field    is   greater    than   that    for 

m 

the  YT  field.   The  thermal  field  Y_  is  zero  in  the  initial  state  but 

it  changes  as  a  result  of  the  Ekman  friction  term.   Y   is  introduced 

m 

into  the  model  as  a  random  disturbance  in  y.   Introduction  of  surface 

friction  leads  to  the  coupling  of  the  two  fields.   Thus  at  time  t  =  0 

A  =  random  disturbance  (y) ,  (5.3) 

B  =  random  disturbance  (y) ,  (5.4) 

C  =  0  (5.5) 

D  =  0  (5.6) 

E  =  Uoyoln(cosh(y/yo)),  (5.7) 

or  E  =  U  y  tanh  y/y  .  (5.8) 

o  o        o 
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6.   COMPUTATIONAL  PROCEDURE 

A  list  of  the  constants  used  in  the  computations  is  given  below; 

y  =  200  kilometers, 
o 

w  =  2000  or  4000  kilometers, 

U  =  10  meters  per  second, 
o  r 

f  =  5  x  10    per  second, 
o  r 

G     =0.8  meters    per    second   squared    per  centibar   squared, 

8      =   0.0   or   2.29   x    10  per  meter    per   second, 

Ap  =  50  centibars, 

A  =  10  meters  per  second  squared, 

a  =  22 .5  degrees , 

2  -12 

|j.  =  2.49  x  10    per  meter  squared, 

H  =  25  or  50  kilometers, 

At  =  0.5  or  1.0  hours . 

When  w  =  2000  km  and  H  =  50  km  a  forecast  period  of  14  days  was  used, 

in  all  other  cases  the  forecast  period  was  28  days.   A  sample  program 

as  developed  by  Schminke  (1968)  and  modified  to  fit  the  IBM  360  computer 

is  included  as  Appendix  1  (note:   this  program  is  for  w  =  2000  km  and 

H  =  25  km). 

The  random  disturbance  introduced  into  the  Y   field,  (5.3)  and 

m 

(5.4),   was    generated   by  a    random   number  generation   program,    then   read 
directly   into   the    computer.      The    same   disturbance    field  was    used    for 
all   experiments  with   both   variations    of   6,    H,    At,    and   all  variations 
of   the    non-dimensional  wave    number.      After   a    few  days    the   disturbance 
adjusted    to  an  exponential  growth.      For   convenience,    growth    rates   were 
computed   over   the    last    three   days    of   the    forecast    period.      Once    the 
field   has   adjusted    it   does    not   matter  when   computed,  rates   are   always 
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the  same.   This  was  done  by  assuming  the  amplitudes  could  be  written  as 

xL  =  Sexp  (UQn/yo)  t^  (6.1) 

x2  =  Sexp  (UQn/yo)  t^  (6.2) 

where  n  is  the  growth  rate  and  S  is  a  constant.  By  forming  the  ratio 

x„/x  ,  the  growth  rate  can  be  shown  to  be 

y   In  (x  /x  ) 

O  Z     1  /£   ON 
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7 .   RESULTS 

A  series  of  experiments  was  conducted  under  varied  conditions1.   In 
general  the  total  distance  over  the  grid  was  2000  km,  but  one  series  was 
done  over  a  total  distance  of  4000  km.   Grid  mesh  was  either  50  km  or  25 
km  with  the  smaller  distance  producing  better  results.   When  the  grid 
mesh  was  25  km  a  time  step  of  0.5  hours  was  used,  however  the  majority 

of  experiments  were  run  at  t  =  1.0  hours.   All  cases  were  run  with 

-11 
3=0,  then  with  p  =  2.29  x  10    meters  per  second.   Many  combinations 

of  friction,  wave  number  and  wind  profile  were  run,  and  selected  ones 
are  illustrated  in  figures  4  through  19. 

Growth  rates  for  the  disturbances  with  both  profiles  and  for  various 
non-dimensional  wave  numbers  are  shown  in  figures  4  through  11.   Betchov 
and  Criminale  (1967)  found  a  maximum  growth  rate  of  approximately  0.19 
for  a  non-dimensional  wave  number  of  0.45  using  the  hyperbolic  tangent 
profile.   With  rigid  boundaries,  no  friction  and  a  fairly  large  grid 
size  (H  =  25  km)  through  the  maximum  shear  zone  (width  approximately 
400  km),  the  result  (fig.  7)  was  a  growth  rate  of  0.165  at  a  wave 
number  of  0.45.   According  to  Yanai  and  Nitta  (1968a)  the  best  results 
should  be  obtained  with  a  minimum  of  20  grid  steps  across  the  zone  of 
maximum  shear.   Therefore  the  results  for  the  hyperbolic  tangent  pro- 
file are  considered  good,  and  the  trend  upward  toward  the  Betchov  and 
Criminale  (1967)  results  is  evident  as  the  grid  size  is  decreased 
(compare  figures  5  and  7).   The  inclusion  of  the  p  parameter  had  no 
significant  effect  on  the  growth  rates  giving  only  a  slight  decrease 
of  0.01  or  less. 
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The  hyperbolic  secant  squared  profile  results  were  equally 
encouraging.   Just  as  Betchov  and  Criminale  (1967)  lead  one  to  expect, 
a  relatively  broad  zone  of  maximum  growth  is  found  reaching  from  0.85 
to  0.95  non-dimensional  wave  numbers.   It  is  difficult  to  ascertain  just 
what  maximum  rate  Betchov  and  Criminale  (1967)  found,  however,  results 
of  these  experiments  show  rates  that  are  significantly  smaller  than  the 
rates  for  the  hyperbolic  tangent  profile.   The  inevitable  conclusion  is 
that  the  hyperbolic  secant  squared  profile  is  not  quite  as  unstable  as 
the  hyperbolic  tangent  profile.   As  before,  decreasing  the  grid  size 
in  the  maximum  shear  area  led  to  larger  growth  rates .   Figures  9  and  10 
illustrate  this  very  well. 

In  figures  12  through  15  the  amplitude  and  phase  relationships  of 

Y  ,  Ym  and  a)  are  depicted  for  the  wave  number  of  maximum  growth, 
m   T 

Amplitudes  of  Y  ,  Y_  and  a)  show  maxima  in  the  region  of  greatest 
horizontal  shear.   The  influence  of  the  frictional  term  is  clearly 
seen  in  the  relative  amplitudes  of  Y,  and  Y- .   The  three  fields  show 
close  phase  relationships  and  all  phases  change  rapidly  through  the 
shear  zone.   It  appears  to  be  a  safe  assumption  that  the  temperature 
and  pressure  fields  are  in  close  relation  for  both  levels.   The  omega 
amplitude  (fig.  14)  shows  a  distinct  minimum  in  the  center  of  the 
field  but  sharp  maximum  peaks  on  each  side  in  the  maximum  shear  zone. 
Phase  relations  between  Yo  and  a)  might  almost  be  termed  ideal  with  the 
sinking  motion  north  of  the  shear  zone  about  70  degrees  upwind  of  the 
amplitude  ridge,  placing  the  convergence  in  the  area  where  it  is 
normally  observed  in  easterly  waves. 

Addition  of  the  0  parameter  made  small  insignificant  changes  in  the 
amplitude;  phase  relationships  maintained  the  same  basic  pattern  as 
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without  3.   Variation  in  the  non-dimensional  wave  number  over  the  range 
0.25  to  0.60  left  the  pattern  basically  unchanged  while  producing  the 
growth  rate  changes  discussed  above. 

Figures  16  through  19  illustrate  the  behavior  of  the  hyperbolic 
secant  squared  profile.   Note  that  Y-,  far  exceeds  the  other  components 
in  amplitude  (fig.  16)  and  that  the  maximum  is  symmetrical  around  the 
center  of  the  maximum  shear  zone.   Y~  is  very  small  indicating  that 
friction  is  probably  more  effective  with  the  hyperbolic  secant  squared 
profile.   The  lower  field,  Y  ,  tends  to  lead  the  upper  field  in  west- 
ward movement.   This  appears  reasonable  and  is  observed  in  tropical 
disturbances  of  this  nature  (Riehl,  1954,  pp.  215).   This  gives  the 
entire  system  a  tilt  toward  the  east  with  increasing  height;  this  tilt 
is  not  noticeable  with  the  hyperbolic  tangent  profile. 

The  omega  amplitude  is  very  sharp  peaked  and  symmetric  about  the 
center.   Velocities,  from  about  +  y   to  the  boundaries,  are  very  small. 
In  this  area  the  convergence  appears  to  be  about  100  degrees  out  of 
phase  with  the  trough.   Near  y  =0,  along  the  Intertropical  Conver- 
gence Zone,  the  convergence  actually  leads  the  trough. 

Use  of  a  finite  value  for  the  |3  parameter  made  no  significant 
difference  in  the  results  with  this  profile. 
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Fig.    4.      Growth    rate    for   hyperbolic    tangent    profile, 
w   ■   2000   km,    H   =   50   km,    with    friction. 


23 


170" 


160— 


130' 


0.30 


0.60 


Fig.  5.   Growth  rate  for  hyperbolic  tangent  profile, 
w  =  2000  km,  H  =  50  km,  no  friction. 
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Fig.  6.   Growth  rate  for  hyperbolic  tangent  profile, 
w  =  2000  km,  H  =  25  km,  with  friction 
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Fig.  7.   Growth  rate  for  hyperbolic  tangent  profile, 
w  =  2000  km,  H  =  25  km,  with  no  friction. 
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Fig.  8.   Growth  rate  for  hyperbolic  tangent  profile, 
w  ■  4000  km,  H  =  50  km,  with  friction. 
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Growth    rate    for   hyperbolic    secant    squared    profile, 
w   »   2000   km,    H   -   50   km,    with    friction. 
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Fig.  10.   Growth  rate  for  hyperbolic  secant 

squared  profile,  w  ■  2000  km,  H  ■  25  km, 
with  friction. 
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Fig.  11.   Growth  rate  for  hyperbolic  secant 
squared  profile,  w  ■  4000  km, 
H  ■  25  km,  with  friction. 
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Fig.  12.   Amplitudes  of  Y  ,  Y. ,  and  Y-  for  hyperbolic  tangent 

profile  plotted  against  arbitrary  scale,  w  =  4000  km, 
H  ■  50  km,  and  wave  number  ■  0.45. 
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Fig.    13.      Phase   relationships   corresponding   to  amplitudes   of   fig.    12 
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Fig.    14.      Amplitude    of   omega,    arbitrary   scale    for    same    case   as    figs 
12  and    13. 
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Fig.  15.   Phase  relationship  between  Y.  and  o>,  same  conditions 
as  figs.  12,  13,  and  14. 
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Fig.  16.   Amplitude  of  Y.,  Y~  ,  and  Y  with  hyperbolic 

secant  squared  profile,  w  ■  2000  km,  H  =  25  km, 

ky  ■  0.95,  and  scale  arbitrary. 

J  o 


34 


Phase   relationship  for   Y.  ,    Y- ,    Y_  for  same  conditions 
as    fig.    16. 
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Fig.  18.   Amplitude  of  omega,  conditions  as  in  fig.  16, 
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Fig.  19.   Phase  relationship  between  omega  and  Y~  for  conditions 
in  fig.  16. 
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8.   CONCLUSIONS 

The  easterly  wave  model  developed  appears  to  be  reasonable  if  the 
assumption  of  barotropic  instability  as  the  energy  source  is  accepted. 
Basic  zonal  flow  has  been  introduced  by  two  profiles,  both  independent 
of  height.   The  temperature  variations  arose  from  the  action  of  the 
Ekman  friction  terms.   Waves  developed  by  the  model  behaved  in  many 
respects  as  easterly  waves  have  been  observed  to  act.   Results  show 
that  the  growth  rate  of  disturbances  is  a  function  of  wave  number,  there' 
fore  wave  length,  as  well  as  other  factors.   Separation  of  the 
boundaries  to  greater  distances  made  some  small  changes  but  the 
greater  effect  is  obtained  by  decreasing  the  grid  size,  which  also 
required  changing  the  time  step  size,  through  the  shear  zone. 

Vertical  motions  produced  by  the  hyperbolic  tangent  profile  are 
very  similar  to  those  observed  in  easterly  waves,  in  that  the  maximum 
sinking  motion  is  upwind  of  the  lower  level  ridge  and  the  maximum 
rising  motion  is  upwind  of  the  lower  level  trough.   The  disturbance 
maintains  its  phase  relationship  in  the  vertical  so  that  no  tilt  is 
observed.   With  the  hyperbolic  secant  squared  profile  the  amplitude 
of  the  vertical  motion  outside  the  maximum  shear  zone  is  small  and 
again  the  rising  motion  is  upwind  of  the  trough  and  sinking  upwind 
of  the  ridge.   Within  the  area  of  large  horizontal  shear  the  con- 
vergent area  moves  to  a  position  downwind  of  the  low  level  trough. 
It  is  not  known  whether  this  is  observationa lly  verified  along  the 
Intertropical  Convergence  Zone.   Using  this  profile,  the  hyperbolic 
secant  squared,  the  model  has  developed  a  solution  that  indicates  a 
definite  eastward  tilt  with  height.   This  tilt  has  been  verified  by 
observations.   Perhaps  some  combination  of  the  two  profiles  would  be 

more  realistic  in  depicting  actual  waves. 
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A  more  complete  investigation  of  the  effect  of  friction  on  the 
symmetric  secant  squared  profile  appears  warranted.   There  is,  also, 
a  definite  need  to  include  the  effects  of  condensation  in  the  model. 

Non-linear  effects,  their  role  in  changing  the  mean  wind  profile 
and  the  structure  of  the  disturbance,  are  very  important  and  must  be 
considered.   Use  of  the  primitive  equations  or  the  balance  equations 
would  seem  to  offer  the  most  effective  means  of  incorporating  these 
effects.   They  appear  to  be  crucial  if  an  easterly  wave  is  to  grow 
into  a  tropical  storm. 

Always  of  utmost  importance  is  the  need  for  more  detailed  and 
extensive  observations,  but  nowhere  is  this  as  true  as  of  the  tropics 
No  model  can  be  developed  and  adequately  tested  without  detailed 
observations . 
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